Directed percolation describes lifetime and growth of turbulent puffs and slugs 
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We show that Directed Percolation (DP) simulations in a pipe geometry in 3+1 dimensions fully 
capture the observed complex phenomenology of the transition to turbulence. At low Reynolds 
numbers (Re), turbulent puffs form and spontaneously relaminarize. At high Re, turbulent slugs 
expand uniformly into the laminar regions. In a spatiotemporally intermittent state between these 
two regimes of Re, puffs split and turbulent regions exhibit laminar patches. DP also captures 
the main quantitative features of the transition, with a superexponentially diverging characteristic 
lifetime below the transition. Above the percolation threshold, active (turbulent) clusters expand 
into the inactive (laminar) phase with a well-defined velocity whose scaling with control parameter 
(Reynolds number or percolation probability) is consistent with experimental results. Our results 
provide strong evidence in favor of a conjecture of Pomeau. 
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The transition from laminar to turbulent flow, studied 
first by Reynolds [1 in pipes, remains a source of complex 
and fascinating phenomenology. Varying the dimension- 
less parameter bearing his name. Re = ULjv^ where U 
and L are characteristic flow velocity and length scales 
and V is the kinematic viscosity, Reynolds observed local- 
ized clusters of turbulence, now commonly called "puffs" , 
that can spontaneously split or decay. Later, Wygnan- 
ski et al [2] systematically described the phase diagram 
of the laminar to turbulent transition as a function of the 
Reynolds number, and the dynamics of the transition is 
an active area of investigation today [3]-[5] . Laminar pipe 
flow is known to be linearly stable for all Reynolds num- 
bers, but small disturbances can trigger a transition to 
the turbulent state [6]. For sufficiently low Re the fluid 
flow is always laminar and any turbulent disturbances 
decay immediately. However, when Re ~ 2000 and the 
inlet disturbance is large enough, metastable turbulent 
puffs can appear. These are regions of turbulent flow of 
a fixed size (on the order of pipe diameter D) that can 
spontaneously relaminarize. 

The lifetime of metastable turbulent puffs has been 
carefully measured for 1650 < Re < 2050, and found 
to grow superexponentially with Re, spanning eight or- 
ders of magnitude in time [7 . For larger values of Re, 
the characteristic lifetime of these puffs grows, and they 
begin to split and show complex spatiotemporal behav- 
ior [8[[9]. In this regime, the puffs interact with a char- 
acteristic distance: if puffs come too close to each other, 
they may spontaneously annihilate [10 . The splitting 
process continues, until Re exceeds a critical value Rcc ~ 
2500, above which and for a sufficiently large inlet distur- 
bance, a uniform state of turbulence, a "slug" , grows with 
clearly defined turbulent-laminar interface and a velocity 
that scales approximately with v^Re — Rcc [JX JJJ • The 
phase diagram of pipe flow turbulence is shown schemat- 
ically in Fig. [l] 

The purpose of this Letter is to show in detail that the 



phenomenology and quantitative details of many features 
of the laminar-turbulence transition are consistent with 
the non-equilibrium phase transition in the universality 
class of directed percolation (DP), as originally conjec- 
tured by Pomeau [13 and continued in subsequent works 
by many authors (see, e.g. [14 and references therein). 
Our work measures the lifetime of active states in DP 
in a pipe geometry, finding agreement with the super- 
exponential functional dependence recently measured by 
Hof et al [7] . We also measure the growth rate of active 
DP clusters in the supercritical directed percolation and 
show that our scaling results are in good agreement with 
available experimental data on the growth rate of turbu- 
lent slugs [11] . These results show that dynamical phase 
transition phenomena may be described by directed per- 
colation, supporting earlier detailed observations of the 
DP critical exponents in a fluctuating turbulent liquid 
crystal system driven by external forcing [151 116]. 

The analogy between DP and turbulent-to-laminar 
transition is the following: Active states in a three di- 
mensional lattice correspond to coarse-grained regions of 
size ~ T] where the turbulence intensity exceeds a thresh- 
old (here 77 is the viscous scale), whereas inactive states 
correspond to patches of the fluid which are laminar. The 
dimension along the percolating direction is associated 
with time t in the usual interpretation of DP as a dynam- 
ical process. The percolating probability p is analogous 
to Reynolds number Re in the vicinity of the percolation 
transition, but the mapping need not be linear. For the 
metastable puffs. Re < 2050 region is mapped to p < Pc 
whereas for the growing fronts. Re > 2500 is mapped to 
p > Pc The critical region maps into the spatiotempo- 
ral regime, as summarized in Fig. [T] but this region and 
Pc is not strictly defined except in the limit of infinite 
system size. We simulate DP in 3+1 dimensions, where 
the inlet disturbance in pipe flow experiments is modeled 
as the initial region of active sites. The simulations are 
performed in the reference frame of the traveling puff. 



2 



metastable 
Re puffs 



splitting puffs/ 
spatiotemporal 
intermittency 



expanding 
slugs 




metastable 
clusters 



criticality 



expanding 
clusters 



FIG. 1. (Color online). Comparison of the phenomenology of 
transitional turbulence as a function of Re with that of DP in 
3+1 dimensions as a function of p, both in a pipe geometry. 




FIG. 2. (Color online) (a) Illustration of the decay of an active 
cluster in subcritical directed percolation in 3+1 dimensions 
in pipe geometry. In this figure the length of the pipe extends 
horizontally, and the active sites in percolation are displayed 
with green cubes. Inactive sites are not shown, (b) Different 
front shapes in the growing front bond DP model. Rough 
fronts occur for small p — Pc whereas smoother fronts occur 
for large p — pc- (c) Time evolution of an initially active 180 
contiguous sites percolating downwards simulated via bond 
percolation for p = 0.61 < Pc- The active state (in black) 
fully decays into the absorbing state (white) after about 250 
steps. 



which usually travels slower than the laminar mean flow 
velocity U. 

Lifetime of turbulent puffs: The survival probability 
of turbulent puffs in pipe flow is known to be memory- 
less [HI HE], 



P(Re, t) = exp 



t-tp 
"r(Re) 



(1) 



where t > to. Here, the survival probability P(Re,t) 
refers to the probability that the turbulent puff still ex- 



ists after flowing for time t, to is the formation time of 
the puff and r(Re) is the characteristic lifetime. In Hof et 
aVs work to was a constant (70 D/U where U is the mean 
flow velocity and D is the pipe diameter) 0. By mea- 
suring P(Re, t) for specific times t and Reynolds num- 
bers Re, Hof et al. calculated r(Re) from ([T]). They dis- 
covered that r scales superexponentially with Reynolds 
number [7 , fitting well a parameterization of the form: 



r(Re) = To exp[exp(ciRe + C2)], 



(2) 



where tq ~ D/U. 

The survival probability P(p, t) in our DP measure- 
ments is the probability that there will be active sites 
left in the lattice after t DP steps. From P{p,t)^ the 
lifetime of the disturbance r can be measured just as in 
Hof et al. This idea is illustrated in the snapshots of the 
simulation in Fig. [2|a). Here, DP is simulated in 3+1 
dimensions in a pipe of radius of 5 lattice sites. In this 
simulation p is less than Pc, and so the puff eventually de- 
cays. One can measure the lifetime with a 3-dimensional 
lattice where two of the spatial dimensions span a disk 
of radius of R lattice points (corresponding to the pipe 
radius), with fixed boundary conditions. However, the 
measurement of lifetime in this way over many orders of 
magnitude is made difficult because of the system size. 
The lifetime measurements must be repeated many times 
to be able to obtain sufficient statistics. However, when 
p is close to the critical percolation threshold Pc^ the cor- 
relation length ^± ~ ip — PcY^ along a space dimension 
becomes larger than i?, and the nominally 3+1 dimen- 
sional DP is effectively 1+1 dimensional. Thus, to get 
sufficient statistics we simulate DP in a 1-dimensional 
lattice of length N that is initially made to be active in 
a subregion of length A^o- Below the DP critical point 
the active states will eventually decay into the absorbing 
state. In a finite-sized system, the decay can always oc- 
cur, but we find that the characteristic lifetime r of the 
decay grows super-exponentially with percolation thresh- 
old p, and beyond a certain percolation probability, the 
average lifetime of the active state is too large to be mea- 
surable on a computer. 

Hof et al. were able to calculate r via ([T]) by measur- 
ing P(Re,t). Even though they could only extend t to 
3450I)//7, they were able to resolve P to 100 ppm, giving 
them effective measurements of r over 8 orders of mag- 
nitude. In the case of directed percolation, one cannot 
use this procedure, since to is not constant, but instead 
depends on the percolation probability p. In directed 
percolation, to is the time over which the initial state is 
remembered by the system. Hence, for each value p we 
must measure the survival lifetimes of many instantia- 
tions of directed percolation. The cumulative distribu- 
tion function (CDF) of these survival lifetimes then ap- 
proximates the survival probability P(p, t), as long as suf- 
ficiently many instantiations have been performed. From 
the fit of the form ([T]) to the CDF data, one can read off 
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FIG. 4. (Color online) Measured values of front propagation 
velocity G (red crosses) compared to theoretical prediction 
for 3+1 dimensional and 1+1 dimensional DP. Green vertical 
line indicates the value of p for which exceeds \.2R. We 
call this value ^r. (a) Plot indicating the power law crossover 
of G(j)). Inset shows that Gjij) — Pc)°"^^^ is roughly constant 
for p < PR and similarly G/ (p — Pc)^'^^^ for p > pR. (b) Both 
regimes of DP have the same critical point pc (within the error 
of our measurements). The linear fits to G'^^^ are shown in 
solid lines. Extrapolations are indicated with dashed lines, 
and they cross at the same pc. This value of pc is corresponds 
to that of 3+1 dimensional bond percolation. 



FIG. 3. (Color online) (a) Superexponential scaling of the 
characteristic lifetime r. The line indicates the fit to Er- 
ror bars indicate 95% confidence intervals from Kolmogorov- 
Smirnov test (for p > 0.62) and 90% confidence intervals from 
test (for p < 0.62). In the inset, tq = 0.017. (b) Numerical 
data for survival probabilities P(p, t) (points) and a fit to ([T]) 
with T given by Q (solid lines). Data shown for t = 300 
(red crosses), t = 1000 (green dots), t = 6000 (blue plusses) 
and t = 80000 (cyan triangles). The value of pc observed cor- 
responds to that of 1+1 dimensional bond percolation, (c) 
Measured survival probabilities as functions of t for 4 differ- 
ent values of percolation probability p. Blue line indicates 
measured data, whereas red line indicates a fit to the expo- 
nential distribution. Deviations from exponential distribution 
for small t are due to nonzero to- 

r{p) and to{p). Our measurements of r{p) for a lattice 
of size N = 100 and Nq = 20 are given in Fig. IsFa). The 
line in the Figure is obtained by fitting r to (l2|). The 
inset shows the linear fit to loglogr(p)/ro. Sufficiently 
far away from the critical point we find that the linear 
fit deviates, indicating that the superexponential behav- 
ior may somehow be related to the diverging correlation 
lengths at the critical point. 

From the CDF data, we can evaluate the survival prob- 
ability functions analogously to Figure 2 of |7j. Fig- 
ure |3jb) shows our numerical data P(p, t) for 4 differ- 
ent times (S curves), along with the model P{p^t) — 
exp( — (t — t^) jrij))). To evaluate the model fit, we use 
the value of ^0(^0.2) where po.2 is found by finding where 
P{Po.2^t) = 0.2. Note that the fact that the S curves 
become steeper with p is a characteristic of the superex- 
ponential scaling of r{p). 

The numerical data presented so far in the paper has 
been measured in a finite volume of size N = 100. Finite 
size effects in DP have been investigated thoroughly in 



the literature [19]. We ran our simulations in a volume 
bound only by the range of the integers on our computer 
(0 to 2^^ — 1), and we didn't find any qualitative differ- 
ences regarding the superexponential scaling of r{p). 

Growth rate model: When p > Pc active DP clusters 
grow in the pipe. We measured their growth rate and 
related it to the growth rate of turbulent slugs. The 
speed at which the front of the percolating clusters prop- 
agates into the neighboring inactive region is given by 
G - ^x/^ii - b-Pc)^i'"^^, where ^ - {p - Pc^ is the 
correlation length in direction of space (denoted by +) 
or in direction of time (denoted by ||). Using the above 
prescription and numerical values of DP critical expo- 
nents [19 one should expect that G ^ {p — Pc)^ where 
7 = 0.524 in 3+1 dimensional DP, whereas 7 = 0.637 
in 1+1 dimensions. These power laws are close to the 
exponent 0.5 first proposed in 1986 [11^, as well as in 
modern experiments [12]. However, the data is not suf- 
ficient yet to differentiate between two such close power 
law exponents. 

Measurements of the growth rate G of an initially ac- 
tive region in 3+1 DP in a pipe geometry are shown 
in Fig. [4j The measurements were made by simulating 
bond DP with p > Pc and measuring the positions of 
the two fronts as functions of time. During the numeri- 
cal simulation we also measure the correlation length 
by calculating the root-mean-square height (i.e. rough- 
ness) of the turbulent-laminar interface. The agreement 
with theoretical expectation is good, and we see numer- 
ical evidence for the crossover from 3+1 to 1+1 dimen- 
sions. When p — pc <^ I, then (,_\_ > R and the system is 
effectively 1+1 dimensional. In Fig. [4]^b), we have plot- 
ted G^/^ versus p for different choices of 7 correspond- 
ing to 1+1 and 3+1 dimensional DP. We see clearly the 
crossover between the expected regimes. Note that in this 



4 



plot we did not need to guess Pc'. both scaling regimes 
yield the same Pc- It is difficult however to extend the 
data for G{p) close to Pc- Due to the finite size of the 
system, when p — pc is small, the active regions split and 
may decay into the absorbing state. This makes it dif- 
ficult to clearly measure front propagation velocity. On 
the other hand, when p — Pc is large, the scaling breaks 
down. Thus we expect the power law exponent of 0.524 
to be observable only in an intermediate regime of p — pc^ 
sufficiently close to Pc but still such that £^ < R. 

One other aspect of the phenomenology of pipe flow is 
captured by the DP model, namely that the fronts of ac- 
tive regions with p — 1 are much rougher than when 
p—Pc is large. This is because the density of active states 
within the region is an increasing function of p. Fur- 
thermore, the width of the front is related to the spatial 
correlation length ^± which becomes small when p — Pc 
is large. The difference between the rough and smooth 
front regimes is shown in Fig. |2|b) . This is analogous to 
the results in pipe flow experiments hot wire measure- 
ments, where puff structures were found to have rough 
edges whereas slugs have clearly defined fronts |2J [20] . 

Conclusion: The DP simulations of characteristic life- 
time presented in this paper have been performed via 
the bond percolation algorithm. However, we found su- 
perexponential scaling of r{p) for site percolation too. 
Therefore, the superexponential scaling of the lifetime 
is likely a universal characteristic of the directed per- 
colation process. Goldenfeld et al proposed that this 
superexponential character of the turbulent puff lifetime 
can be described by extreme value statistics [2r , because 
puff decay occurs when turbulent energy fails to attain 
the required threshold at all points in the puff [22]. In 
the usual central limit theorem, under appropriate condi- 
tions [23 the distribution of a sum of N random variables 
tends to a Gaussian distribution for large N. However, 
a maximum (or minimum) of N random variables is in- 
stead superexponentially distributed with 3 universality 
classes [21] , selected by the underlying probability distri- 
bution of {xi}. In Fig. [2|c) we show a time evolution of 
an initially active cluster percolating with p < pc. The 
lifetime of the entire cluster is the lifetime of the longest 
active "strand" percolating downwards. Assuming that 
strand lifetimes are independent and identically (expo- 
nentially) distributed, then the lifetime of the longest 
strand is given by the type I Fisher-Tippett distribution 
exp(— exp(— p)). This argument has also been used to ex- 
plain the superexponential distribution of size of largest 
connected cluster in ordinary (isotropic) percolation [24 , 
and it was found there that correlations between cluster 
sizes (analogous to strand lifetimes) do not influence the 
superexponential scaling. 

The DP model we proposed in this paper can account 
for the superexponential lifetime of the turbulent puffs, 
as well as the uniform growth rate of turbulent slugs. 
As shown in Fig. [T] the transition between these two 



regimes (2050 < Re < 2500) occurs through the splitting 
and interactions of puffs. The spatio-temporal patterns 
of coarse-grained turbulent intensity obtained from a di- 
rect numerical simulation ^ bear similarities to those 
of directed percolation, but the data are not adequate 
to make a quantitative analysis. Thus, we do not know 
if additional hydrodynamic realism is required for a de- 
tailed account of the regime 2040 < Re < 2500. After 
this work was complete, we learned of a new preprint by 
Barkley, which uses a coupled map and a shear field to 
model the laminar-turbulent transition [25] . It is possible 
that this model is in the DP universality class also. 
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